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Abstract 

Motivated by the very recent measurement of the branching ratio of B® — > J/ipg decay, we 
calculate the branching ratios of Bj — ► J/ipn and B^ — > J/^rj' decays in the perturbative 
QCD (pQCD) approach. The pQCD predictions for the branching ratios of considered decays 
are: BR(B® — > J/^/g) = (1.96lg'g|) x 10~ 6 , which is consistent with the first experimental 
measurement within errors; while BR(B% -> J/Vrf) = (1.09±g;||) X 10 , very similar with 
— ► J/^g decay and can be tested by the forthcoming LHC experiments. The measurements 
of these decay channels may help us to understand the QCD dynamics in the corresponding 
energy scale, especially the reliability of pQCD approach to these kinds of B meson decays. 
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Very recently, the first observation of B® — > J/^/t] decay was reported by Belle Collab- 
oration jlj, and the branching ratio measured is 

BR(B° d -> J/Vrj) = (9.5 ± 1.7 (stat) ± 0.8(syst)) x 1(T 6 , (1) 

which is consistent with the currently available theoretical predictions [H, @, 0] . 

Up to now, the theoretical calculations for the branching ratios of B^ — > J/^rj^ decays 
were obtained by using the heavy quark factorization approximation in Ref. 0], or from 
the measured J/^/7r° and J/^K° branching ratiosjl, 0,1^] based on the assumption of 
the SU(3) flavor symmetry of strong interaction. In this paper, we will calculate the 
branching ratios of B® — > J/tyrj and B® — > J/^rf'^ decays directly by employing the 
low energy effective Hamiltonian 0] and the perturbative QCD (pQCD) factorization 
approach @, H, 0| • 

The paper is organized as follows: we present the formalism used in the calculation of 
B® — > J/ipr]^ decays in Sec. [I] In Sec. [TTJ we show the numerical results and compare 
them with the measured values. A short summery and some conclusions are also included 
in this section. 



I. FORMALISM AND PERTURBATIVE CALCULATIONS 

The pQCD approach has been developed earlier from the QCD hard-scattering ap- 
proach ML and has been used frequently to calculate various B meson decay channels 
@, 0, 0, LL3]. F° r two body charmless hadronic B^ s —> Mrft (here M stands for the 
pseudo-scalar or vector light mesons composed of the light quarks u,djs) decays, the 
pQCD predictions generally agree well with the measured values 



In Refs. [H, (3, the authors calculated B -> D*K, D { * )+ D ( * ] and B s -> £>(*>+£>(*)- 
decays and found that the pQCD approach works well for such decays. Here we try to 
apply the pQCD approach to calculate the B meson decays involving the heavier J/ty 
meson as one of the two final state mesons. 



A. Formulism 

In pQCD approach, the decay amplitude of B — > J/tyP (P = r),r)('> here) decay can 
bo written conceptually as the convolution, 

A(B -> MxM 2 ) ~ Jd%d%d% Tr [Cit^B^^j/^k^p^H^^k^t)] , (2) 

where the term "Tr" denotes the trace over Dirac and color indices. C(t) is the Wilson 
coefficient which results from the radiative corrections at short distance. In the above 
convolution, C(t) includes the harder dynamics at larger scale than Mb scale and describes 
the evolution of local 4-Fermi operators from (the W boson mass) down to t ~ 
0(a/AMb) scale, where A = Mb — rrib. The function H(ki, k2, k%, t) is the hard part and 
can be calculated perturbatively. The function $m is the wave function which describes 
hadronization of the quark and anti-quark to the meson M. While the function H depends 
on the process considered, the wave function <3>m is independent of the specific process. 
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Using the wave functions determined from other well measured processes, one can make 
quantitative predictions here. 

Using the light-cone coordinates the B meson and the two final state meson momenta 
can be written as 

Mr , . M B , 9 . Mn , o . 

P x = -^(1,1, Or), P 2 = -|(l,r 2 ,0 T ), P 3 = -|(0,l-r 2 ,O r ), (3) 

respectively, where r = Mj/q, /Mb, and the light meson masses mjj have been ne- 
glected. The longitudinal polarization vector of the J/^f meson, ex, is given by = 
^=^f — (1, —r 2 ,Or)- Putting the light (anti-) quark momenta in B, J/^> and rf) mesons 

as k\, k 2 , and fc 3 , respectively, we can choose 

h = ( Xl P+,0,k 1T ), k 2 = (x 2 P+,0,k 2T ), k 3 = (0,x 3 P 3 -,k 3T ). (4) 

Then, for B — > J/^frj decay for example, the integration over k^ , k^ , and & 3 in eq. (121) 
will lead to 



A{B -> J/Vrf) ~ J i 



dxidx2dx^bidbib2db2b^db^ 

•Tr [C(t)$ B (x u 6i)$j/*(a; 2 , & 2 )^(x 3 , & 3 )# &i, ^tfo) e~ 5 W] ,(5) 

where 6j is the conjugate space coordinate of kix, and t is the largest energy scale in 
function H(xi, b i: t). The large logarithms ln(mw/i) are included in the Wilson coefficients 
C(t). The large double logarithms (ln 2 Xj) on the longitudinal direction are summed by 



the threshold resummation [14|, and they lead to S t (xi) which smears the end-point 
singularities on x\. The last term, e -5 *-*-', is the Sudakov form factor which suppresses the 
soft dynamics effectively [15| . Thus it makes the perturbative calculation of the hard 
part H applicable at intermediate scale, i.e., Mb scale. We will calculate analytically the 
function H(xi, bi,t) for the considered decays in the first order in a s expansion and give 
the convoluted amplitudes in next section. 

B. The B% -> J/^lrf) Decays 

The low energy effective Hamiltonian for decay modes P° — > J/iprj^ can be written as 

n eff = 2| [v cb v; d (c^oftO + c 2 ^)o^))} , (6) 

with the four-fermion operators 

0\ = J«7 M (1 - 75)0,3 • c/ry^l - 7 5 )& Q , O c 2 = d a ^(l - 7 5 )c Q • cp^l - 75)6,3 (7) 

where the Wilson coefficients Cj(/x) (i = 1,2), we will use the leading order (LO) expres- 
sions, although the next-to-leading order (NLO) results already exist in the literature 0. 
This is the consistent way to cancel the explicit /1 dependence in the theoretical formulae. 
For the renormalization group evolution of the Wilson coefficients from higher scale to 
lower scale, we use the formulae as given in Ref.[l6[ directly. 
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FIG. 1: Typical Feynman diagrams contributing to the Cabibbo- and color- suppressed B® 
J/tyrjO decays. 



As for B meson wavefunction, we make use of the same parameterizations as used 
in the studies of different processes For vector J /if} meson, in terms of the nota- 

tion in Ref. [j~7j ]. we decompose the nonlocal matrix elements for the longitudinally and 
transversely polarized J /if} mesons into 



m J ,4 L ^ L {x)+/ L I/^ t {x / 



Here, \1/ L denote for the twist-2 distribution amplitudes, and for the twist-3 distri- 
bution amplitudes, x represents the momentum fraction of the charm quark inside the 
charmonium. 



The J /ifj meson asymptotic distribution amplitudes read as [18 

x(l — x) 



V L (x) 



10.94 



2V2N C 

fj/ip 



0.7 



2 v / 2iVc 



2x) [ 



1 - 2.8x(l - x) 
x(l — x) 
1 - 2.8x(l -x) 



0.7 



(9) 



It is easy to see that both the twist-2 and twist-3 DAs vanish at the end points due to 
the factor [x(l — x)] ' 7 . 

From the effective Hamiltonian ([6]), the Feynman diagrams corresponding to the con- 
sidered decay are shown in Fig.l. With the meson wave functions and Sudakov factors, 
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the hard amplitude is given as 

F ev = SnCpm^ / dx\dx 3 / bxdbxb 3 db 3 4>b{ x x, &i) 
Jo Jo 

x {(1 - r 2 ) [(1 + x 3 (l - r 2 ))rf(x 3 , 63) + r (l - 2x 3 ) 

•<(x 3 , 63)] + r [(1 - 20:3) + r 2 (l + 2x 3 )] <(x 3 , 63) 

•o; s (tg) h e (x!, x 3 , 61, 63) expf-S 1 ^^)] 

+2r [1 - (1 - xi)r 2 0^(a; 3 , 63) - x 1 r 2 ^(a;3, 63)] 

■a s (t 2 e )h e (x 3 , x 1 ,b 3 ,b 1 ) exp[-S ab (t 2 e )]} . (10) 

where r = mj/me; Cp = 4/3 is a color factor. The function h e , the scales t\ and the 
Sudakov factors S a b are displayed in Appendix [A] 

For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 63 can be performed using 5 function £(6 3 — 61), leaving only 
integration of 61 and 62. For the concerned operators, the corresponding decay amplitude 
is 

16v^6 f°° 

M erj = hCfto^ I dx\dx 2 dx 3 I b\db\b 2 db 2 0/3 s (^i, 61) 

3 J J Q 

x {2rr c 4>j/^(x 2 , 6 2 )0^(x 3 , 6 2 ) - 4rr r c (f) t j / ^(x 2 , b 2 )4%(x 3 , 6 2 ) 
- [x 2 r 2 + x 3 (l- 2r 2 )] </>5 /9 (x 2 , h)rf(x 3 , 6 2 ) 
+2 [x 3 r + (x 2 - x 3 )r r 2 ] (frj /q ,(x 2 , b 2 )(j)^(x 3 , 6 2 ) 

■a s (t f )h f (x 1 , x 2 , x 3 , 61, 6 2 ) exp[-S cd (t f )]} . (11) 

where r c = m c /mB,m c is the mass for c quark. 

For the B% — > J/^rf decay, the Feynman diagrams are obtained by replacing the 77 
meson in Fig. 1 with the meson rj . The corresponding expressions of decay amplitudes will 
be similar with those as given in Eqs. fllOtiTTI) . since the 77 and 7/ are all light pseudoscalar 
mesons and have the similar wave functions. The expressions of — > J/^/t]' decay can 
be obtained simply by the following replacements 

< — <— > < — r ^r' . (12) 

For the 77-/7' system, there exist two popular mixing basis: the octet-singlet basis and 
the quark- flavor basis [ill, 20]. Here we use the quark- flavor basis Ji| and define 

rj q = (uu + dd)/V2, Vs = ss. (13) 
The physical states r\ and 77' are related to r] q and r] s through a single mixing angle 0, 

V \ _ TT( M\ ( Vq \ _ ( C0S( P -Sm0\ / T] q 



U(4>) [ ,q ) = [ ."I ™y '« . (14) 
V J \Vs J y sin cos <p J \ T] S J 

The three input parameters f g , /. and in the quark-flavor basis have been extracted 
from various related experiments 19, 2^ 

f q = (1.07 ±0.02)/,, /, = (1.34 ±0.06)/,, = 39.3° ± 1.0°, (15) 
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where f n = 130 MeV. In the numerical calculations, we will use these mixing parameters 
as inputs. It worth of mentioning that the effects of possible gluonic component of rf 
meson will not considered here since it is small in size [lOj, |2l|, |22 . 



For B® — > J/^/rj decay, by combining the contributions from different diagrams, the 



total decay amplitude can be written as 



M(B° d -> J/Vri) = VoV^FM Fejj,* 



Ci + -C 2 



M e „C2 



(16) 



where the relevant mixing parameter is F 1 ((f)) = cos0/\/2. 

It should be mentioned that the Wilson coefficients Ci = Ci(t) in Eq. (fTBT) should 
be calculated at the appropriate scale t using equations as given in the Appendices of 
Ref. [l6j. Here the scale t in the Wilson coefficients should be taken as the same scale 
appeared in the expressions of decay amplitudes in Eqs. ffTOj) and fTlTl) . This is the way 
in pQCD approach to eliminate the scale dependence. In order to estimate the effect of 
higher order corrections, however, we introduce a scale factor at — 1.0 ± 0.2 and vary the 
scale t max as described in Appendix A. 

Similarly, the decay amplitudes for B® — > J/^/r/' decay can be obtained easily from 
Eq. ffTBl) by the following replacements of Fi((fi) — > F[((p) = sin <p/y/2. 



II. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will calculate the branching ratios for those considered decay modes. 
The input parameters and the wave functions to be used are given in Appendix [B] In 
numerical calculations, central values of input parameters will be used implicitly unless 
otherwise stated. 

With the complete decay amplitudes, we can obtain the decay width for the considered 
decays, 

Y{B^J/^)^ L ^{1-t^)\m{B^J/^)\ . (17) 

By employing the quark-flavor scheme of 77 — 77' system and using the mixing parameters 
as given in Eq. (TTSJ), one finds the branching ratios for the considered two decays with 
error bars as follows: 



Br( B° d 
Br( B° d -> J/^rj 



J/^rj) = [l.96lS:S(^)lgS(a*)to:SS(a2)lo:S(/j/»)] >< iO" 6 , (18) 
JM) = [l.09^;^(^)^; 7 ?(a 4 ):^?(a 2 )^(/ J/1& )] x 10~ 6 , (19) 

where the main errors are induced by the uncertainties of uo^ = 0.40 ± 0.05 GeV, a t = 
1.0 ± 0.2, a 2 = 0.115 ± 0.115 and f J/9 = 0.405 ± 0.014 GeV , respectively. One can see 
that the pQCD predictions are sensitive to the variations of uj^ and a t . 

For B® — > J /^T] decay, the central value of the pQCD prediction for Br(B d — > J/^Jr]) 
is a factor of 4 smaller than the measured value as given in Eq. ([1]) [lj. But the pQCD 
prediction is in fact still consistent with Belle's first measurement if we take the large 
theoretical and experimental errors into account. By varying the scale factor a t in the 
range of a t = [0.8, 1.0], for example, the central value of Br(B —* J/^/rj) will change in the 
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range of [0.2, 1.1] x 10~ 5 accordingly. It is not difficult to understand such a t dependence. 
Since the J/^f meson is much heavier than light mesons, and therefore moving not as fast 
as those light meson when B meson is decaying. So a small decrease of the scale will 
lead to a larger Wilson coefficients C\^{t) and ol s {U)i and consequently results in a larger 
decay rate. 

For B® — > J/^/rj' decay, only experimental upper limit (at 90% C.L) is available now: 
BR(B° -> J/*t/) < 6.3 x 10" 5 0, Q. The pQCD prediction for the branching ratio of 
B® — > J/tyrf decay is very similar in magnitude with that of — > J/^r], consistent with 
the upper limit and will be tested in the forthcoming LHC experiments. 

At the leading order, only the tree Feynman diagrams as shown in Fig. 1 contribute to 
B® — ► J/i&rf'' decays. There exists no CP violation in these decays within the standard 
model, since there is only one kind of Cabibbo-Kabayashi-Muskawa (CKM) phase involved 
in the corresponding decay amplitudes, as can be seen from eq. f fT6|) . 

In short, we calculated the branching ratios of B® — ► J/^rj and B® — > J/^/rj' decays at 
the leading order by using the pQCD factorization approach. Besides the usual factoriz- 
able diagrams, the non-factorizable spectator diagrams are also calculated analytically in 
the pQCD approach. By keeping the transverse momentum kx, the end-point singularity 
disappears in our calculation. 

From our calculations and phenomenological analysis, we found the following results: 

• Using the quark-flavor scheme, the pQCD predictions for the branching ratios are 



where the various errors as specified previously have been added in quadrature. 

• The major theoretical errors of the pQCD predictions are induced by the uncertain- 
ties of the hard energy scale tiS and the parameters Uj,. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function h^s, coming from the Fourier transformations of the function 




(20) 
(21) 



H(°\ 



h e {x!,x 3 , 6i,6 3 ) = K ( a/xix 3 (1 -r 2 )m B bi) 9(h - b 3 )K ( ^3(1 - r 2 )m B bi 



■I (^3(1-^)77^63) + 0(63 - h)K (y X3 (l - r 2 )m B 6 3 ) 




(Al) 
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h f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = I 6{b 2 - b 1 )l (M B y / x 1 x 3 (l - r 2 )6 1 )K (Af BV /xix 3 (l - r 2 )b 2 ) 



+ (h 



K (M B F (1) b 2 ) 



for F 2 1} > 



f hS\M b J\F*~\ b 2 ), for F 2 } < 



(A2) 



where Jo is the Bessel function, K and Jo are the modified Bessel functions with 
K (—ix) = — (n/2)Y (x) + i(iT/2)J (x), and F^s are defined by 



F (i) = ( x i ~ ^2)2:3(1 - r<2 ) + 
(xi - x 2 )x 3 (l - r 2 ) + r 



F {2) 



c 1 
2 

c ' 



The threshold resummation form factor St(xi) is adopted from Ref. [17 



2l+2cr(3/2 + c) 

t{) ~ 0FT(l + c) [X{1 X)J ' 

where the parameter c = 0.3. This function is normalized to unity. 
The Sudakov factors used in the text are defined as 



(A3) 
(A4) 



(A5) 



S ab (t) 



Scd(t) 



s (xim B / V2, bij + s (x 3 m B / V%, & 3 ) + s ((1 
1 



x 3 )m B 



/V2,b, 



In Ht/A \ I In ln(t/A) 



-ln^A) -ln(6 3 A) 
s (xiin B /V2,b^j +s (x 2 m B /V2,b 2 

+s (x 3 m B / \/2, bij + s ((1 - x 3 )m B / y/2, bij 
1 



(A6) 



s ((1 - x 2 )m B /V2, 6 2 ) 



2 In ln(t/A) +ln ln(t/A) 
-ln(6iA) -ln(6 2 A)_ 



(A7) 



where the function s(g, b) are defined in the Appendix A of Ref. [16j. The scale V s i n the 
above equations are chosen as 



t\ = a t ■ max(v/x 3 (l - r 2 )M B , l/b u l/6 3 ), 
t 2 = at ■ max( v /x 1 (l - r 2 )M B , I/61, l/6 3 ), 

t/ = at • max( v /x 1 x 3 (l - r 2 )M B , y/( Xl - x 2 )x 3 (l - r 2 ) + r 2 M B , I/61, 1/62), (A8) 



where at = 1.0 ±0.2 and r = Mj/^/M B . The scale tj's are chosen as the maximum energy 
scale appearing in each diagram to kill the large logarithmic radiative corrections. 

APPENDIX B: INPUT PARAMETERS AND WAVE FUNCTIONS 



The masses, decay constants, QCD scale and B® meson lifetime are 
A^r 4) = 250MeV, f n = 130MeV, f J/9 = 405MeV, 



rn. 



1.08GeV, M B o = 5.28MeV, M 



J/* 



3.097GeV, 



M w = 80.41GeV, r B o = 1.54 x 10~ 12 s 



(Bl) 
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For the CKM matrix elements, here we adopt the Wolfenstein parametrization for the 
CKM matrix, and take A = 0.2272, A = 0.818, p = 0.221 and 77 = 0.340 0]. 
For the B meson wave function, we adopt the model 



4>B(x,b) = Nbx 2 (1 — x) 2 exp 



M 2 B x 2 
2ul 



:b2) 



where u;& is a free parameter and we take ujb = 0.40 ± 0.05 GeV in numerical calculations, 
and Nb = 91.745 is the normalization factor for Ub = 0.40 for the B meson. 
The wave function for dd components of rf^ meson is given by 



lib 

V2N C 



Vdd K 



(B3) 



where p and x are the momentum and the momentum fraction of r] d g respectively, while 
(p^ dd , - and (j)^ d _ represent the axial vector, pseudoscalar and tensor components of the 
wave function respectively. We here assume that the wave function of ri dd - is same as the 
7r wave function based on SU(3) flavor symmetry. The parameter ( is either +1 or — 1 
depending on the assignment of the momentum fraction x. 



The explicit expression of chiral enhancement scale m\ 
,2 1 



m 



Vdd 



is given by [21j 



q 



m 



2m n 



x r 2 2 
\m„ cos ■ 

2m„ r > 



I 2 • 2 / V^/s / 2 

+ m„/Sin <p — (m„/ 



m v ) cos (psinc, 



and numerically m\ = 1.07MeV for m v = 547.5 MeV, = 957.8 MeV, f q 
f s = 1.34/,, and = 39.3°. 



(B4) 
1.07/,, 



For the distribution amplitude 



£ and (fi^ q , we utilize the results for n meson 



obtained from the light-cone sum rule 23] including twist-3 contributions: 



X 



[5(1 - 2a;) 2 - 1] 



(B5) 



+a?^ [21(1 - 2xf - 14(1 - 2x) 2 + l] } , 

~ ^O^ - \p 2 n )j [3(1 - 2x) 2 - 1] 

\pl q{s) al^j [35(1 - 2x) 4 - 30(1 - 2x) 2 + 3] \ ,(B6) 



2^/2N c 
1 



, 27 2 

3r]3Uj3 ~ 2(A ~ 10 Pr ' 



V2N ( 
1 



=/,(!- 2x) 



( , + (5773 - ifftws - ^ - jjp 2 ^' 8 )(10* 2 - 10s + 1) 



(B7) 



with the updated Gegenbauer moments [24 



0.115, 



-0.015, p Vq = 2m g /m gg , 773 = 0.015, w 3 = -3.0. (B8) 
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